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Abstract 

We present in detail a dispersion relation formalism for virtual Compton scat- 
tering (VCS) off the proton from threshold into the A(1232)-resonance region. 
Such a formalism can be used as a tool to extract the generalized polarizabil- 
ities of the proton from both unpolarized and polarized VCS observables over 
a larger energy range. We present calculations for existing and forthcoming 
VCS experiments and demonstrate that the VCS observables in the energy 
region between pion production threshold and the A(1232)-resonance show 
an enhanced sensitivity to the generalized polarizabilities. 



PACS numbers : 11.55.Fv, 13.40.-f, 13.60.Fz, 14.20.Dh 
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Typeset using REVT^ 



I. INTRODUCTION 



The field of virtual Compton scattering (VCS) lias been opened up experimentally in re- 
cent years by the new high precision electron accelerator facilities. On the theoretical side, 
an important activity has emerged over the last years around the VCS process in different 
kinematical regimes (see e.g. |l|J^ for reviews). 

In VCS off a nucleon target, a virtual photon interacts with the nucleon and a real photon 
is emitted in the process. At low energy of the outgoing real photon, the VCS reaction 
amounts to a generalization of real Compton scattering (RCS) in which both energy and 
momentum of the virtual photon can be varied independently, which allows us to extract 
response functions, parametrized by the so-called generalized polarizabilities (CPs) of the 
nucleon On the other side, VCS has also a close relation to elastic electron scattering. 
More precisely this means, that the physics addressed with VCS is the same as if one would 
perform an elastic electron scattering experiment on a target placed between the plates of a 
capacitor or between the poles of a magnet. In this way one studies the spatial distributions 
of the polarization densities of the target, by means of the CPs, which are functions of the 
square of the four-momentum, Q^, transferred by the electron. The GPs teach us about the 
interplay between nucleon-core excitations and pion-cloud effects, and their measurement 
provides therefore a new test of our understanding of the nucleon structure. 

A first dedicated VCS experiment was performed at the MAMI accelerator, and two com- 
binations of the proton GPs have been measured Q]. Further experimental programs are 
underway at the intermediate energy electron accelerators ( JLab , MIT-Bates , MAMI 
0) to measure the VCS observables. 

At present, VCS experiments at low outgoing photon energies are analyzed in terms of 
low-energy expansions (LEXs). In the LEX, only the leading term (in the energy of the 
real photon) of the response to the quasi-constant electromagnetic field, due to the internal 
structure of the system, is taken into account. This leading term depends linearly on the 
GPs. As the sensitivity of the VCS cross sections to the GPs grows with the photon energy, 
it is advantageous to go to higher photon energies, provided one can keep the theoretical 
uncertainties under control when approaching and crossing the pion threshold. The situa- 
tion can be compared to RCS, for which one uses a dispersion relation formalism to 
extract the polarizabilities at energies above pion threshold, with generally larger effects on 
the observables. 

It is the aim of the present paper to present in detail such a dispersion relation (DR) 
formalism for VCS on a proton target, which can be used as a tool to extract the GPs from 



VCS observables over a larger energy range, into the A(1232)-resonance region. In Ref. [jlO 



we have given a first account of the DR predictions for the GPs. In this paper we present 
the formalism in detail and show the results for the VCS observables. 

In Sec. we start by specifying the kinematics and the invariant amplitudes of the VCS 
process. 

In Sec. |T|, we set up the DR formalism for the VCS invariant amplitudes and show that 
for 10 of the 12 VCS invariant amplitudes unsubtracted DRs hold. 

In Sec. it is shown that the DR formalism provides predictions for 4 of the 6 GPs of 
the proton. 

In Sec. 0, it is discussed how the s-channel dispersion integrals, which correspond to the 
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excitation of vrA^, ttitN,... intermediate states, are calculated. In the numerical evaluation 
of the dispersion integrals, only the contribution of vrA^ states are taken into account. 

In Sec. [VT] , we show how to deal with the two VCS invariant amplitudes for which one 
cannot write down an unsubtracted DR. Our DR formalism involves two free parameters, 
being directly related to two GPs, and which are to be extracted from a fit to experiment. 

In Sec. [VI I| , we show the results in the DR formalism for both unpolarized and polar- 
ized VCS observables below and above pion threshold. We compare with existing data and 
present predictions for planned and forthcoming experiments. 

Finally, we present our conclusions in Sec. IVTTl. 



Several technical details on VCS invariant amplitudes and hehcity amplitudes are collected 
in three appendices. 



II. KINEMATICS AND INVARIANT AMPLITUDES FOR VCS 

In this section, we start by briefly recalling how the VCS process on the proton is accessed 
through the ep — *• epy reaction. In this process, the final photon can be emitted either by the 
proton, which is referred to as the fully virtual Compton scattering (FVCS) process, or by 
the lepton, which is referred to as the Bethe-Heitler (BH) process. This is shown graphically 
in Fig. ^ leading to the amplitude T*^^ of the ep eprf reaction as the coherent sum of 
the BH and the FVCS process : 

rpee'-f _ rpBH _|_ rpFVCS ^-^^ 

The BH amplitude T^^ is exactly calculable from QED if one knows the nucleon electro- 
magnetic form factors. The FVCS amplitude T^^^^ contains, in the one-photon exchange 
approximation, the VCS subprocess 7*p ■yp. We refer to Ref. [Q] where the explicit expres- 
sion of the BH amplitude is given, and where the construction of the FVCS amplitude from 
the 7*p — >• 7p process is discussed. In this paper, we present the details how to construct 
the amplitude for the 7*p — > jp VCS subprocess, in a DR formalism. 

We characterize the four- vectors of the virtual (real) photon in the VCS process 7*p —>■ jp 
by Q {q') respectively, and the four-momenta of initial (final) nucleons by p {p') respectively. 
In the VCS process, the initial photon is spacelike and we denote its virtuality in the usual 
way by = - Q^. Besides Q^, the VCS process can be described by the Mandelstam 
invariants 

s = iq + p)\ t={q-q')\ u = {q - p')\ (2) 

with the constraint 

s + t + u = 2M^ - , (3) 

where M denotes the nucleon mass. We furthermore introduce the variable u, which changes 
sign under s ^ u crossing : 

(4) 



4M 



and which can be expressed in terms of the virtual photon energy in the lab frame {E^ ) as 
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In the following, we choose Q^, v and t as the independent variables to describe the VCS 
process. In Fig. ^ we show the Mandelstam plane for the VCS process at a fixed value of 
= 0.33 GeV^, at which the experiment of |^ was performed. 
The VCS helicity amplitudes can be written as 

TvAW;AA. = -e%(g,A) e':{q',X') il{p' , X'^) M^'' u{p, X^) , (6) 

with e the proton electric charge (e^/47r = 1/137.036). The polarization four-vectors of the 
virtual (real) photons are denoted by e {e ), and their helicities by A (A'), with A = 0, ±1 
and A' = ±1. The nucleon helicities are Xn,X'^ = ±1/2, and u,u are the nucleon spinors 
(as specified in appendix |Cp. The VCS tensor J^^^'^ in Eq. (^ can be decomposed into a 
Born (B) and a non-Born part (NB) : 

A^/^- = + M%'s . (7) 

In the Born process, the virtual photon is absorbed on a nucleon and the intermediate 
state remains a nucleon, whereas the non-Born process contains all nucleon excitations and 
meson-loop contributions. The separation between Born and non-Born parts is performed 
in the same way as described in Ref. [Q], to which we refer for details. 
One can proceed by parametrizing the VCS tensor of Eq. (|^) in terms of 12 independent 



amplitudes. In Ref. [|1T|, a tensor basis was found so that the resulting non-Born invariant 
amplitudes are free of kinematical singularities and constraints, which is an important prop- 
erty when setting up a dispersion relation formalism. In detail, we denote the tensor Ai^^" 
as Q 



12 



-M^"- = E MQ^'^,t)pf\ (8) 



i=l 



where the 12 independent tensors pf*^ are given in appendix ^ The 12 independent invariant 
amplitudes /« are expressed in terms of the invariants Q^, v and t, but are otherwise identical 
with the amplitudes used in |Tl|. 
The tensor basis p^" of Eq. (X2) was chosen in [11| such that the resulting invariant 



amplitudes /, are either even or odd under crossing. Photon crossing leads to the symmetry 
relations among the /j at the real photon point : 

/,(Q2 = o,i.,t) = +/,(Q2 = 0,-J.,t) , (2 = 1,2,6,11), 

/,(Q2 = o,i/,t) = -/,(Q2 = 0,-z/,t) , (2 = 4,7,9,10), (9) 

while the amplitudes /s, /s, /§, /12 do not contribute at the real photon point, because the 
corresponding tensors in Eq. (|A2|) vanish in the limit —>■ 0. 

Nucleon crossing combined with charge conjugation provides the following constraints on 
the fi at arbitrary virtuality : 

h (g^ z/, t) = + /, {Q\ -V, t) , {t = 1, 2, 5, 6, 7, 9, 11, 12) , 

/, (Q^ z/, t) = - /, {Q\ -z/, t) , (z = 3, 4, 8, 10) . (10) 
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When using dispersion relations, it will be convenient to work with 12 amplitudes that are 
all even in v. Therefore, we define new amplitudes Fi {i = 1,...,12) as follows : 



1,2,5,6,7,9,11,12) 
= 3,4,8,10), 



(11) 



satisfying F^ (Q^, — z/, t) = F^ (Q^, z^, t) for i = 1,...,12. As the non-Born invariant amplitudes 
/34^8 10 ^ ^ z/ — i> 0, the definition of Eq. ( ]TT| ) ensures that also all the non-Born Fl'^ {i 
= 1,...,12) are free from kinematical singularities. The results for the Born amplitudes F^^ 
are listed in Appendix ^. 

From Eqs. (|^) and ([T0|), one furthermore sees that F^ and Fg vanish at the real photon 
point. Since 4 of the tensors vanish in the limit 0, only the six amplitudes Fi, F2, F4, 

Fq, Fio and Fn enter in real Compton scattering (RCS). 

Dispersion relation formalisms for RCS were worked out in Refs. pjj^ in terms of another 
set of invariant amplitudes, also free from kinematical singularities and constraints and de- 
noted as Ai{h', t) [i = 1,...,6) (see Appendix A of Ref. p for definitions). It is therefore useful 
to relate the amplitudes -Fi,2,4,6,io,ii(0, i^, t) to the RCS amplitudes (z/, t) (? = 1, 6). We 
find after some algebra the following relations at = : 



-e^ F2 



-e^F, 



-e'F, 



-A, - 
1 

' A 



A. 



4M2 



A, 
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2M2 
1 

4M2 
1 

~2M 
1 

~4M 



4 , 



t-4M^ 



+ 



4M2 

[^5 - A^\ , 

t - 4M2 + 4z/2 



A^ 



4M2 



Aa + A. 



(12) 



where the charge factor — appears explicitely on the Ihs of Eq. ([T^), because this factor is 
included in the usual definition of the A^. The values of the RCS invariant amplitudes Ai {% 
= 1,...,6) at z/ = t = can be expressed in terms of the scalar polarizabilities a, /3, and the 
spin polarizabilities 71,72,73,74, as specified in Ref. f^. 



III. DISPERSION RELATIONS AT FIXED T AND FIXED fqR VCS 

With the choice of the tensor basis of Eq. (|A^), and taking account of the crossing relation 
Eq. (|13), the resulting non-Born VCS invariant amplitudes Fi {% = 1,...,12) are free of all kine- 
matical singularities and constraints and are all even in z/, i.e. -Fi(Q^, z^, t) = Fi{Q'^, —i^,t). 

Assuming further analyticity and an appropriate high-energy behavior, the amplitudes 
Fi{Q'^, z/, t) fulfill unsubtracted dispersion relations with respect to the variable u at fixed t 
and fixed virtuality : 
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ReFr{Q\ iy,t) = - V du' 7^"^ ' ^ , (13) 
where we indicated explicitely that the Ihs of Eq. (|13D represents the non-Born {NB) parts of 



the amphtudes. Furthermore, in Eq. ([T3|) , Im^Fj are the discontinuities across the s-channel 
cuts of the VCS process, starting at the pion production threshold, which is the first inelastic 
channel, i.e. uthr = f^n + ijni + + Q^/2)/(2M), with m^r the pion mass. 

Besides the absorptive singularities due to physical intermediate states which contribute 
to the rhs of dispersion integrals as Eq. (0), one might wonder if other singularities exist 
giving rise to imaginary parts. Such additional singularities could come from so-called 
anomalous thresholds 0,0] , which arise when a hadron is a loosely bound system of other 



hadronic constituents which can go on-shell (such as is the case of a nucleus in terms of 
its nucleon constituents), leading to so-called triangular singularities. It was shown that in 
the case of strong confinement within QCD, the quark-gluon structure of hadrons does not 
give rise to additional anomalous thresholds [^,0, and the quark singularities are turned 
into hadron singularities described through an effective field theory. Therefore, the only 
anomalous thresholds arise for those hadrons which are loosely bound composite systems of 
other hadrons (such as e.g. the S particle in terms of A and vr). For the nucleon case, such 
anomalous thresholds are absent, and the imaginary parts entering the dispersion integrals 
as in Eq. (0) are calculated from absorptive singularities (due to vrA^, vrvrA^, ... physical 
intermediate states). 

The assumption that unsubtracted dispersion relations as in Eq. ( p!3|) hold, requires that 
at high energies {v ^ oo aX fixed t and fixed Q^) the amplitudes \m.sFi{Q'^, u, t) {i = 1,...,12) 
drop fast enough so that the integrals of Eq. (jl^) are convergent and the contribution from 
the semi-circle at infinity can be neglected. 

For the RCS invariant amplitudes Ai,...,Aq which appear on the rhs of Eq. (p!^, Regge 
theory leads to the following high-energy behavior for u oo and fixed t : 

Ai,A2 ~ , (14) 

{As + A,) ~ z/"-W-2 ^ (15) 

^3,^5 ~ zy"M(*)-2 (16) 
~ Zy"M(t)-3 ^ (17) 

where ^^/(t) < 0.5 (for t < 0) is a meson Regge trajectory, and where ap{t) is the Pomeron 
trajectory which has an intercept ap(0) ~ 1.08. Note that the Pomeron dominates the high 
energy behavior of the combination of A^ + Aq. From the asymptotic behavior of Eqs. (|l^ - 
T7D , it follows that for RCS unsubtracted dispersion relations do not exist for the amplitudes 



Ai and A2. The reason for the divergence of the unsubtracted integrals is essentially given 
by fixed poles in the t-channel, notably the exchange of the neutral pion (for A2) and of 
a somewhat fictitious a- meson (for Ai) with a mass of about 600 MeV and a large width, 
which models the two-pion continuum with the quantum numbers I = J = 0. 

We consider next the VCS amplitudes Fi, ...,Fi2, in the Regge limit {u ^ 00 at fixed t 
and fixed Q"^) to determine for which of the amplitudes unsubtracted dispersion relations 
as in Eq. (0) exist. The high-energy behavior of the amplitudes Fi is deduced from the 
high-energy behavior of the VCS helicity amplitudes that are defined and calculated in 
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Appendix |C[ This leads, after some algebra, to the following behavior in the Regge limit 
(z/ — ^ oo, at fixed t and fixed Q"^) Q : 







^ap{t)-2 




(18) 


Fs + 4Fn 




^ap{t)-2 




(19) 


F2, Fg, Fio 




jj»p{t)-2 ^ 


jy<^M(t)-2 ^ 


(20) 


F7 








(21) 


F3, Fg 






jjaM{t)~2 


(22) 


Fg, F12 




^ap{t)-i ^ 


jjaM{t)-2 ^ 


(23) 


F4 




^ap{t)-i 




(24) 



In Eqs. (0 - 1^), we have indicated the high energy behavior from the Pomeron (ap) and 
from the meson (a^) contributions separately. It then follows that for the two amplitudes 
Fi and Fs, an unsubtracted dispersion integral as in Eq. (|13D does not exist, whereas the 
other ten amplitudes on the Ihs of Eqs. (^ - ^) can be evaluated through unsubtracted 
dispersion integrals as in Eq. (p!3|). 

Having specified the VCS invariant amplitudes and their high energy behavior, we are now 
ready to set up the DR formalism. First, we will show in Sec. |IV| that 4 of the 6 GPs of the 
nucleon can be evaluated using unsubtracted DR. We will then discuss in Sec. how the 
s-channel dispersion integrals of Eq. (|TB|) are evaluated. In particular, unitarity will allow us 
to express the imaginary parts of the VCS amplitudes in terms of ttN, ttttN,... intermediate 
states. Finally, we will show in Sec. [V^ how to deal with the remaining two VCS invariant 
amplitudes for which one cannot write unsubtracted DRs. 



IV. DISPERSION RELATIONS FOR THE GENERALIZED POLARIZABILITIES 

The behavior of the non-Born VCS tensor M.'^^ '^^ -^l- (H) energy (q' = |g'| — > 0) 

but at arbitrary three- momentum q = |g| of the virtual photon, can be parametrized 
by six generalized polarizabilities (GPs), which are functions of q and which are denoted 
by pip' L',pL)S(^^^^ p,p^,pi|. In this notation, p (p') refers to the electric (2), magnetic (1) 
or longitudinal (0) nature of the initial (final) photon, L [L' = 1) represents the angular 
momentum of the initial (final) photon, and S differentiates between the spin-flip [S = 1) 
and non spin-flip {S = 0) character of the transition at the nucleon side. A convenient choice 
for the 6 GPs has been proposed in : 

p(oi,oi)0(q)^ P(ii'ii)0(q), (25) 
p(oi,oi)i(q)^ p(ii,ii)i(q)^ piW^^^^ P(0i,i2)i^q^_ (26) 



^We note that some of the Fi in Eqs. (IS - ^) decrease faster with increasing v than reported 
in Ref. [I0|. This is because a more detailed calculation has shown a cancellation in the highest 
power of u for some of the Fj, which leads to the behavior of Eqs. ( [l8| - 24). However, this does 
not change the conclusion obtained in Ref. [^] that unsubtracted DR only exist for 10 of the 12 
Fi. The asymptotic behavior of Eqs. (18 - 24) only shows that for some of those 10 amplitudes, 
the dispersion integrals converge even faster than anticipated earlier |1C|. 
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In the limit q — for the GPs, one finds the following relations with the polarizabilities 



(in gaussian units) of RCS [[TT 



p(01,0l)0(^Q^) 



4:71 2 



An 

^ V 3 

4:71 y/2 



p(oi-oi)i(o) = o, 

p(ii,ii)i(o) = o, (27) 

In terms of invariants, the limit q' — > at finite three-momentum q of the virtual photon 
corresponds to — *• and t —Q^ at finite . One can therefore express the GPs in terms 
of the VCS invariant amplitudes Fj at the point = 0,t = —Q^ for finite Q^, for which we 
introduce the shorthand : 

F.(Q2) ^ F^^[Q\v = ^,t = -Q^) . (28) 

The relations between the GPs and the FiiQ"^) can be found in . 



The present work aims at evaluating the GPs through unsubtracted DRs of the type of 
Eq. ([T3|). We have seen from the high-energy behavior that the unsubtracted DRs do not 
exist for the amplitudes Fi and F5, but can be written down for the other amplitudes. 
Therefore, unsubtracted DRs for the GPs will hold for those GPs which do not depend 
on the two amplitudes Fi and F^. However, the amplitude F5 can appear in the form 
F5 -f 4Fii, because this combination has a high-energy behavior (Eq. (pO])) leading to a 
convergent integral. Among the six GPs we find four combinations which do not depend on 
Fi and F5 : 

+ V"'")° = ^ (^^±^y^'Mgo {|f2 + (2F6 + F9)-Fi2}, (29) 

P^°^'°^)^ = ^ (^^)''%~o{(F5 + F, + 4Fn)+4MF,,}, (30) 

p(oi,i2)i _ 1 p(ii,ii)i ^ 1 ( E + m V" Mqo 
V2qo ^\ E J q2 

X { (F5 + + 4 Fii) + 4 M (2 Fe + i^g) } , (31) 



p(oi,i2)i _^ V^p{ii,02)i ^ i / + ^ \ Qo 



2 ' Q\ E J 

X {go(F5 + F7 + 4Fn) +8M2(2F6 + F9)}, (32) 



where E = + denotes the initial proton cm. energy and qo = M — E the virtual 
photon cm. energy in the limit q' = 0. For small values of q, we observe the relation 
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go ~ ~q^/(2iVf). Furthermore, in the hmit q' = 0, the value of is always understood as 
being = — Q'o, which we denote by for simplicity of the notation. 

The four combinations of GPs on the Ihs of Eqs. (|29| - ^2|) can then be evaluated in 
a framework of unsubtracted DRs through the following integrals for the corresponding 

FAQ'') ■ 



TT /,, W 



V. S-CHANNEL DISPERSION INTEGRALS 



The imaginary parts of the amplitudes Fi in Eq. ([13|) are obtained through the imaginary 
part of the VCS helicity amplitudes defined in Eq. (|^). The latter are determined by using 
unitarity. Denoting the VCS helicity amplitudes by Tfi, the unitarity relation takes the 
generic form 

2 Im, Tf, = Y,{'^^Y^\Px - P■:)T]c^Tx^ , (34) 

X 

where the sum runs over all possible intermediate states X. In this work, we are mainly 
interested in VCS through the A(1232)-resonance region. Therefore, we restrict ourselves 
to the dominant contribution by only taking account of the vrA^ intermediate states. The 
influence of additional channels, like the ttttN intermediate states which are indispensable 
when extending the dispersion formalism to higher energies, will be investigated in a future 
work. 

The VCS helicity amplitudes can be expressed by the Fi in a straightforward manner, 
even though the calculation is cumbersome. The main difficulty, however, is the inversion 
of the relation between the two sets of amplitudes, i.e., to express the twelve amplitudes Fi 
in terms of the twelve independent helicity amplitudes. To solve this problem we proceeded 
in two different ways. First, the inversion was performed numerically by applying differ- 
ent algorithms. Second, we succeeded in obtaining an analytical inversion using a two-step 
procedure. To this end we used an additional set of amplitudes, called Bi, which were intro- 
duced by Berg and Lindner ||T^ and which are defined in Appendix |A 2|. Both the relations 



between the Bi and the Fi on the one hand, and between the helicity amplitudes and the Bi 
on the other hand can be inverted analytically. The expressions of the Fi amplitudes in terms 
of the Bi amplitudes are given in Appendix [A 2[ and the expressions of the Bi amplitudes 



in terms of the VCS helicity amplitudes are given in Appendix |C 3| (for the definition of the 
VCS helicity amplitudes, see Appendices |C 1| and |C 2|) . In our calculations, we checked that 



the two methods to express the Fi amplitudes in terms of the VCS helicity amphtudes lead 
numerically to the same results. 

The imaginary parts of the s-channel VCS helicity amplitudes are calculated through uni- 
tarity taking into account the contribution from vrA^ intermediate states. They are expressed 
in terms of pion photo- and electroproduction multipoles as specified in Appendix |C4[ For 
the calculation of the pion photo- and electroproduction multipoles, we use the phenomeno- 
logical MAID analysis pO| , which contains both resonant and non-resonant pion production 
mechanisms. 
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VI. ASYMPTOTIC PARTS AND DISPERSIVE CONTRIBUTIONS BEYOND ttN 

To evaluate the VCS amplitudes Fi and in an unsubtracted DR framework, we proceed 
as in the case of RCS This amounts to perform the unsubtracted dispersion integrals ( p!3[ ) 
for Fi and F^ along the real z/-axis only in the range —Vmax < ^ < +i^max-, and to close the 
contour by a semi-circle with radius Vmax in the upper half of the complex z/-plane, with the 
result 

ReFf''{Q\ V, t) = Ft\Q^ t) + F^\Q\ z/, t) , (for 2 = 1, 5) , (35) 
where the integral contributions F?"* (for i = 1, 5) are given by 

2 fi'max J.I T™ p. (Q2 I J.\ 

Fr^iQ^ u,t) = du' ^^^^^W,^,^) (f^^ ^ ^ 5) ^ (3g) 

and with the contributions of the semi-circle of radius Vmax identified with the asymptotic 
contributions (-Ff*, F^^). 

Evidently, the separation between asymptotic and integral contributions in Eq. ( p5D is 
specified by the value of Vmax- The total result for Ff'^ is formally independent of the 
specific value of Vmax- In practice, however, v^^x is chosen to be not too large so that 
one can evaluate the dispersive integrals of Eq. ( |5BD from threshold up to v^ax sufficiently 
accurate. On the other hand, Vmax should also be large enough so that one can approximate 
the asymptotic contribution F"* by some energy-independent (i.e. //-independent) function. 
In the calculations, we therefore choose some intermediate value v^^x ~ 1-5 GeV, and 
parametrize the asymptotic contributions F"* by t-channel poles, which will be discussed 
next for the cases of F^^ and F'^^ . 



'as 



A. The asymptotic contribution Fi 

The asymptotic contribution to the amplitude F^ predominantly results from the t-channel 
7r°-exchange, 

FtiQ',u,t) ^ Ff{Q\t) = -AF^l{Q\t) = _L ^-^^^-^-^7(^') . (37) 

As mentioned before, the 7r°-pole only contributes to the amplitudes F5 and Fn, but drops 
out in the combination (F5 + AFn), which therefore has a different high-energy behavior 
as expressed in Eq. (|1^). In Eq. (|37|), the nNN coupling g^NN is taken from Ref. pT 



9^^^/(4:71) = 13.73. Furthermore, in (^Tf), F^o^-y (Q^) represents the 7r''7*7 form factor. Its 
value at = is fixed by the axial anomaly : Ft^o^^ (0) = 1/(4 tt^/^) = 0.274 GeV~^, 
where = 0.0924 GeV is the pion decay constant. For the Q^-dependence of F^o^^ (Q^), 
we use the interpolation formula proposed by Brodsky-Lepage [^] : 

^^"^^ ^^'^ = l + QViSn^P^) ' ^^^^ 



which provides a rather good parametrization of the vr 7*7 form factor data over the whole 

^ 



range, and which leads to the asymptotic prediction at large : F^o^^ {Q^ ^) 
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When fixing the asymptotic contribution F^^ through its vr^-pole contribution as in 
Eq. p7|), one can determine one more GP of the nucleon, in addition to the four com- 
binations of Eqs. - ^2]). In particular, the GP p(^i'^^)^ can be expressed by : 



P^''^''^HQ') = [^^y^^im') + <ioFMQ')}. (39) 

In Fig. ^, we show the results of the dispersive contribution to the four spin GPs, and compare 
them to the results of the 0{p^) heavy-baryon chiral perturbation theory (HBChPT) []23| , 



the linear a-model [Q, and the nonrelativistic constituent quark model |2^. It is obvious 
that the DR calculations show more structure in than the different model calculations. 

The 0{p^) HBChPT results predict for the GPs p(oi'Oi)i and p("'ii)i a rather strong in- 
crease with Q^, which would have to be checked by a O^p"^) calculation. 

The constituent quark model calculation gives negligibly small contributions for the GPs 
p(oi,oi)i g^j^^ p(ii,02)i^ whereas the GPs p(ii'ii)i and p(oi'i2)i receive their dominant contri- 
bution from the excitation of the A(1232) (Ml Ml transition) and N* and A* resonances 
{El M2 transition) respectively. 

The linear a-model, which takes account of part of the higher order terms of a consistent 
chiral expansion, in general results in smaller values for the GPs than the corresponding 
calculations to leading order in HBChPT. 

The comparison in Fig. |^ clearly indicates that a satisfying theoretical description of the 
GPs over a larger range in is a challenging task. 

In Fig. ^, we show the dispersive and 7r°-pole contributions to the 4 spin GPs as well as 
their sum. For the presentation, we multiply in Fig. ^ the GPs p(oi'i2)i ^^d p(ii'02)i -with 
Q, in order to better compare the dependence when including the vr'^-pole contribution, 
which itself drops very fast with Q^. The vr^-pole does not contribute to the GP p(oi.oi)i^ 
but is seen to dominate the other three spin GPs. It is however possible to find, besides the 
GP p(oi'Oi)^^ the two combinations given by Eqs. ( |31|J3^ ) of the remaining three spin GPs, 
for which the vr^-pole contribution drops out [jlO| . 



B. The asymptotic part and dispersive contributions beyond vrA^ to Fi 

We next turn to the high-energy contribution to Pi. As we are mainly interested in a 
description of VCS up to A(1232)-resonance energies, we saturate the dispersion integrals 
by their ttN contribution. Furthermore, we will estimate the remainder by an energy- 
independent function, which parametrizes the asymptotic contribution (i.e. the contour 
with radius i^max in the complex z/-plane), and all dispersive contributions beyond the vrA^ 
channel up to the value Umax = 1-5 GeV. 

Before turning to the case of VCS, we briefly outline the parametrization of the asymptotic 
part of Pi in the case of RCS, and how one expresses it in terms of a polar izability, which 
is then extracted from a fit to experiment. 

The asymptotic contribution to the amplitude Pi originates predominantly from the t- 
channel vrvr intermediate states, and will be calculated explicitly in two model calculations. 
In the phenomenological analysis, this continuum is parametrized through the exchange of a 
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scalar-isoscalar particle in the t-channel, i.e. an effective "(T"-meson, as suggested in Ref. 0. 
For RCS, this leads to the parametrization of the difference of F^^ and its vrA^ contribution, 
as an energy-independent function : 

Fr(Q2 = 0,i.,t)-F-^(Q2 = 0,i.,t) ^ [Fr(0,0,0)-F-^(0,0,0)] ^ _ , (40) 

where F^^ on the Ihs and rhs are evaluated through a dispersive integral as discussed in 
section 0. In Eq. (|40|) , the effective "a" -meson mass uicj is a free parameter in the RCS 
dispersion analysis, which is obtained from a fit to the t-dependence of RCS data, and turns 
out to be around ~ 0.6 GeV 0. The value Ff^{0, 0, 0) is then considered as a remaining 
gobal fit parameter to be extracted from experiment. It can be expressed physically in terms 
of the magnetic polarizability (3 : 



F,^^(0,0,0) = ^P. (41) 



An 



e 



In RCS, one usually takes (a — /?) as fit parameter instead of P because the sum (a + P) 
at the real photon point can be determined independently, and rather accurately, through 
Baldin's sum rule, which leads for the proton to the phenomenological value [^] : 

a + p = (13.69 ± 0.14) x lO^^fml (42) 

Using a dispersive formalism as outlined above, the most recent global fit to RCS data for 
the proton yields the following values for the electric and magnetic polarizabilities of the 
proton : 

a = (12.1 ± 0.3(stat.) T 0.4 (syst.) ± 0.3 (model) ) x 10-^fm^ (43) 
/? = (1.6 ± 0.4(stat.) ± 0.4 (syst.) ± 0.4 (model) ) x 10-^fm^ (44) 

From Eqs. (^, ^), one then obtains for the difference {a — P), the following global average 
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a-P = (10.5 ± 0.9 (stat. + syst.) ± 0.7 (model) ) x 10-^fm^ (45) 

The term Ff^(0,0,0) in Eq. (^0]), when calculated through a dispersion integral, has the 
value : 

Fi-^(0,0,0) = ^p-^ = ^ (9.1 X . (46) 

From the vrA^ contribution P'^^ of Eq. (|i6D, and the phenomenological value P of Eq. (PD, 
one obtains the difference : 

iP - /3"^) = -7.5 X 10"^ fm^ , (47) 

which enters in the rhs of Eq. (140|) . By comparing the value of Eq. (|47D with the total value 
for P (Eq.(P^), one sees that the small experimental value of the magnetic polarizability 
comes about by a near cancellation between a large (positive) paramagnetic contribution 
{P'^'^) and a large (negative) diamagnetic contribution {P — P'^'^), i.e. the asymptotic part 
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of Fi parametrizes the diamagnetism. 

Turning next to VCS, we proceed analogously by parametrizing the non-Born term 
F(^^{Q'^,i>,t) beyond its nN dispersive contribution, by an energy-independent t-channel 
pole of the form : 



^, t) - F-^(Q^ t) ^ , (48) 



where the parameter mo- is taken as for RCS : mo- ~ 0.6 GeV. The function /(Q^) in Eq. ( ^HD 
can be obtained by evaluating the Ihs of Eq. ( ^81) at the point where the GPs are defined, 
i.e. V = {] and t = —Q^, at finite Q^. This leads to : 

fm = [Fim-F^'^m] (i+qv^^) , (49) 

where we introduced the shorthand Fi{Q'^) as defined in Eq. (pH]). Fi{Q'^) can be expressed 
in terms of the generalized magnetic polarizability P'^^^'^^^°((5^) of Eq. ( [25| ) as |Tl|] : 

HQ') = - v/f (^^) P^'^'^nQ') (50) 



4:71 ( IE 



2 \E + M 



1/2 



, (51) 



where 13{Q'^) is the generalized magnetic polarizability, which reduces at = to the 
polarizability of RCS. 
Eqs. (|¥8| , then lead to the following expression for the VCS amplitude 

pNB . 

F^iQ^i^.t) ^ F^iQ^v^t) + [F,{Q^) - F^iQ')] \t.T/"^' ' (52) 

where the vrA^ contributions F^^{Q^,u,t) and F['^{Q'^) (or equivalently /3''^(Q^)) are cal- 
culated through a dispersion integral as outlined above. Consequently, the only unknown 
quantity on the rhs of Eq. (|5^ ) is Fi^Q"^), which can be directly used as a fit parameter 
at finite Q^. This amounts to fit the generalized magnetic polarizability (3{Q'^) from VCS 
observables. 

The parametrization of Eq. (|52D for Fi permits to extract j3{Q'^) from VCS observables at 
some finite and over a larger range of energies with as few model dependence as possible. 
In the following, we consider a convenient parametrization of the dependence of (3{Q'^) 
in order to provide predictions for VCS observables. For this purpose we use a dipole form 
for the difference of /5(Q^) — /^'^^(Q^), which enters in the rhs of Eq. (0) via Eq. (pTI). This 
leads to the form : 



where the RCS value {(3 — (3'"^) on the rhs is given by Eq. (^T]). The mass scale in Eq. 
determines the dependence, and hence gives us the information how the diamagnetism is 
spatially distributed in the nucleon. Using the dipole parametrization of Eq. (|53|) , one can 
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extract from a fit to VCS data at different values. 

To have some educated guess on the physical value of A/3, we next discuss two microscopic 
calculations of the diamagnetic contribution to the GP (3{Q'^). The diamagnetism of the 
nucleon is dominated by the pion cloud surrounding the nucleon. Therefore, we calculate 
the diamagnetic contribution through a dispersion relation estimate of the t-channel tttt 
intermediate state contribution to Fi. Such a dispersive estimate has been performed before 
in the case of RCS where it was shown that the asymptotic part of Fi can be related 



to the 77 — > TTTT — > A^A^ process. The dominant contribution is due to the tttt intermediate 
state with spin and isospin zero (/ = J = 0). The generalization to VCS leads then to the 
identification of F^'^ with the following unsubtracted DR in t at fixed energy z/ = : 

prm - r (54) 

J Ami t' + 



The evaluation of the imaginary part on the rhs of Eq. (^^, originating mainly from the 
TTTT intermediate state contribution, requires information on the subprocesses 7*7 tttt 
and TTTT — i> A^A^. For the latter we use the extrapolation of Ref. |]2B| for the vrA^-scattering 



amplitude to the unphysical region of positive t. For the 7*7 tttt amplitude, we use the 
unitarized Born amplitude, following Ref. 0. At the pion electromagnetic vertex, the pion 
electromagnetic form factor is included. At = 0, it was found that the unitarization 
procedure enhances the 77 tttc cross section in the threshold region, compared to the 
Born result, which is required to get agreement with the data. This becomes obvious from 
the DR of Eq. (|5^) , where the imaginary part of Fi is weighted by 1/t, so that the thresh- 
old contribution dominates the dispersion integral. The dispersive evaluation of Eq. ( |5^ 
contains no free parameters as it uses as input the 77 tttt and tttt NN processes, and 
therefore provides a more microscopic model for the phenomeno logical "a" -exchange. For 
RCS, the dispersion integral Eq. (|^) yields the value ~ —7.3 x 10~^fm^. However, the 
unsubtracted dispersion integral can only be evaluated up to —t = 0.778 GeV^, because the 



TTTT NN amplitudes of Ref. were only determined up to this value, and the disper- 
sion integral of Eq. (^) may not have fully converged at this value. Therefore, one should 
consider the near perfect agreement between the value of from this calculation with the 
phenomenological value of (|47|) as a coincidence. However, our estimate indicates that the 
dispersive estimate through tttt t-channel intermediate states provides the dominant phys- 
ical contribution to the diamagnetism, and that it can be used to give a first guess of the 
distribution of diamagnetism in the nucleon. With this model we show the dependence 
of " in Fig. |. 

To have a second microscopic calculation for comparison, we also show in Fig. |^ an evalua- 
tion of in the linear a-model (LSM) of Ref. |24|. The LSM calculation overestimates 



the value of Ff'*(0) (or equivalently (3as) by about 30% at any realistic value of (which 
is a free parameter in this calculation). However, as for the dispersive calculation, it also 
shows a steep dependence. 

Furthermore, we compare in Fig. ^ the two model calculations discussed above with the 
dipole parametrization for /3(Q^) — [3'^^ {Q'^) of Eq. ( ^3]) for the two values : A^ = 0.4 GeV 
and A^ = 0.6 GeV. It is seen that these values are compatible with the microscopic esti- 
mates discussed before. In particular, the result for A^ = 0.4 GeV is nearly equivalent to 
the dispersive estimate of tttt exchange in the t-channel. The value of the mass scale A^ is 
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small compared to the typical scale of ^ 0.84 GeV appearing in the nucleon magnetic 
(dipole) form factor. This reflects the fact that diamagnetism has its physical origin in the 
pion cloud, i.e. is situated in the surface region of the nucleon. 



C. Dispersive contributions beyond vrA^ to F2 



Though we can write down unsubtracted DRs for all invariant amplitudes (or combinations 
of invariant amplitudes) except for Fi and F5, one might wonder about the quality of our 
approximation to saturate the unsubtracted dispersion integrals by irN intermediate states 
only. We shall show that this question is particularly relevant for the amplitude F2, for 
which we next investigate the size of dispersive contributions beyond the vrA^ channel. We 
start with the case of RCS, where one can quantify the higher dispersive corrections to F2, 
because the value of F<^^ at the real photon point can be expressed exactly (see Eqs. (|2l 
9f)) in terms of the scalar polarizability sum (a + /3) as : 



Fr{0,0,0) = -- 



An 1 
2 (2M)2 



[a + (3) 



The TiN dispersive contribution to (a + /?) provides the value : 

{a + py^ = 11.6 X 10-^fm^ 



(55) 



(56) 



which falls short by about 15 % compared to the sum rule value of Eq. (|42|) . The remaining 
part originates from higher dispersive contributions {mrN, ...) to F2. These higher dispersive 
contributions could be calculated through unitarity, by use of Eq. (|3^) , similarly to the nN 
contribution. However, the present data for the production of those intermediate states (e.g. 
7*A^ ttttN) are still too scarce to evaluate the imaginary parts of the VCS amplitude F2 
directly. Therefore, we estimate the dispersive contributions beyond vrA^ by an energy- 
independent constant, which is fixed to its phenomenological value at u = t = 0. This 
yields : 



FnQ' = o,u,t) 



Fl 



N 



{2My 



TZN 



(57) 



which is an exact relation at z/ = t = 0, the point where the polarizabilities are defined. 

The approximation of Eq. (^) to replace the dispersive contributions beyond vrA^ by a 
constant can only be valid if one stays below the thresholds for those higher contributions. 
Since the next threshold beyond vriV is mrN, the approximation of Eq. ( |57| ) restricts us in 
practice to energies below the A(1232)-resonance. If one wanted to extend the DR formalism 
to energies above two-pion production threshold, one could proceed in an analogous way by 
replacing Eq. (pTf) as follows : 



Fi^^{Q' = 0,u,t) 



tp-kN 

Att 



(g^ = 0,z.,t) + Fr^(Q2 



{2My 



ttN 



(a + P) 



ttttN 



(58) 



i.e. the energy-dependence associated with nN and nnN dispersive contributions would have 
to be calculated explicitly and the remainder be parametrized by an energy-independent 
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constant fixed to the phenomenological value of {a + P). Eq. (0), and Eq. ( ^OD for 

pNB 

modified in an analogous way to include the niiN dispersive contributions, would then allow 
an extension of the DR formalism to energies into the second resonance region. Such an 
extension remains to be investigated in a future work, but because of the present lack of 
experimental input for the ttttN channel, we restrict ourselves in the present work to energies 
up to the A(1232)-resonance region. 

We next consider the extension to VCS, and focus our efforts to describe VCS into the 
A(1232)-resonance region. Analogously to Eq. (|57D for RCS, the dispersive contributions 
beyond vrA^ are approximated by an energy- independent constant. This constant is fixed at 
arbitrary Q^, u = 0, and t = —Q^, which is the point where the GPs are defined. One thus 
obtains for F^^ : 

F2^^(Q^z/,t) ^ F,-^(Q^z/,t) + [F,{Q') - F^W')] , (59) 

where F2{Q'^) is defined as in Eq. (^), and can be expressed in terms of GPs. In this paper, 
we saturate the three combinations of spin GPs of Eqs. - ^) by their ttN contribution, 
and calculate the fourth spin GP of Eq. ( p9| ) through its vrA^ contributions plus the 7r°-pole 
contribution as shown in Fig. ^. Therefore, we only consider dispersive contributions beyond 
the ttN intermediate states for the two scalar GPs, which are then two fit quantities that 
enter our DR formalism for VCS. In this way, and by using Eq. (^9]), one can write the 
difference F2{Q'^) — F|'^((5^) entering in the rhs of Eq. (^) as follows : 

X { [am - a-^m] + [pm - r^(Q')] } , (eo) 

where (3{Q'^) is the generalized magnetic polarizability of Eq. ([51|). Furthermore, ^(Q^) is 
the generalized electric polarizability which reduces at = to the electric polarizability 
a of RCS, and which is related to the GP P^^^'^^^^{Q'^) of Eq. (0) by : 

p(oi,oi)0(Q2) ^ J_L^_^^Q2y (61) 

We stress that Eqs. ( p2D and ( ^9]) are intended to extract the two GPs a(Q^) and l3{Q^) 
from VCS observables minimizing the model dependence as much as possible. As in the 
previous case for /5(Q^), we next consider a convenient parametrization of the dependence 
of a(Q^) in order to provide predictions for VCS observables. Again we propose a dipole 
form for the difference a(Q^) — a'^'^{Q'^) which enters in the rhs of Eq. (pO|), 

where the dependence is governed by the mass-scale A^, again a free parameter. In 
Eq. (PD, the RCS value 



(a-a"^) = 9.6 X 10-^fm^ (63) 

is obtained from the phenomenological value of Eq. (^) for a, and from the calculated vrA^ 
contribution : a'"^ = 2.5 x 10~^fm^. Using the dipole parametrization of ([6^), one can then 
extract the free parameter Aq, from a fit to VCS data at different values. 
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VII. RESULTS FOR EP EPj OBSERVABLES AND DISCUSSION 

Having set up the dispersion formalism for VCS, we now show the predictions for the 
different ep epy observables for energies up to the A(1232)-resonance region. The aim of 
the experiments is to extract the 6 GPs of Eqs. ( p5| , p6D from both unpolarized and polarized 
observables. We will compare the DR results, which take account of the full dependence 
of the ep — > ep7 observables on the energy (q') of the emitted photon, with a low-energy 
expansion (LEX) in q'. In the LEX of observables, only the first three terms of a Taylor 
expansion in q' are taken into account. 

In such an expansion in q', the experimentally extracted VCS unpolarized squared ampli- 
tude A^'^P takes the form || : 

_^exp ^ ^ ^ _^cxp ^ ^^^^ ^g^^ 

q' q' 

Due to the low-energy theorem (LET), the threshold coefficients Al^l" and TW^f are known 
(see Ref. for details). The information on the GPs is contained in A4'^^, which contains 
a part originating from the (BH-|-Born) amplitude and another one which is a linear combi- 
nation of the GPs, with coefficients determined by the kinematics. It was found in Ref. P] 
that the unpolarized observable A4'^'^ can be expressed in terms of 3 structure functions 
Pll{(i), ^Tr(q), and PLriq) by : 

Mr - ^o^^-"^"™ = 2i^2 [ePLLiq) - PTr(q)] + [v^ - ^25 (1 + £)P^r(q)| , (65) 

where K2 is a kinematical factor, e is the virtual photon polarization (in the standard nota- 
tion used in electron scattering), and vi, V2, V3 are kinematical quantities depending on e and 
q as well as on the polar and azimuthal angles (O^^'™' and $, respectively) of the produced 
real photon (for details see Ref. [Q). 

After some algebra, one finds that the 3 unpolarized observables of Eq. ( p3| ) can be ex- 
pressed in terms of the 6 GPs as : 

P^^ = -2v/6MGi=;P(°^'°^)° , (66) 

P^^ = -3Gm-(p^'''''^' - V2qoP^'''''^') , (67) 

Plt = xH^G, P(^i'^^)° + I^Gm P^'''''^' , (68) 
V ^ V ^ qo 

where Ge and Gm stand for the electric and magnetic nucleon form factors Ge{Q^) and 
GuiQ"^), respectively. 

In Fig. |, we show the calculations of Pll — Ptt/^ and Plt, which have been measured at 
MAMI at = 0.33 GeV^ 0. The virtual photon polarization e is fixed to the experimental 
value [e = 0.62), and for the electromagnetic form factors in Eqs. (p6|-|68|) we use the Hohler 
parametrization ||3^ as in the analysis of the MAMI experiment 0]. 

In the lower panel of Fig. ^ the Q^-dependence of the VCS response function Plt is 
displayed, which reduces to the magnetic polarizability j3 at the real photon point (Q^ = 
0). At finite Q^, it contains both the scalar CP (3{Q'^) and the spin CP p(oi'Oi)i^ as seen 
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from Eq. (|68|). It is obvious from Fig. ^ that the structure function Pit results from a 
large dispersive vrA^ contribution and a large asymptotic contribution (to f3) with opposite 
sign, leading to a relatively small net result. At the real photon point, the small value of f3 
is indeed known to result from the near cancellation of a large paramagnetic contribution 
from the A-resonance, and a large diamagnetic contribution (asymptotic part). The latter 
is shown in Fig. ^ with the parametrization of Eq. (|53D for the values A/j = 0.4 and A/j = 
0.6 GeV, which were also displayed in Fig. ^ Due to the large cancellation in Plt, its 
dependence is a very sensitive observable to study the interplay of the two mechanisms. In 
particular, one expects a faster fall-off of the asymptotic contribution with in comparison 
to the ttN dispersive contribution, as discussed before. This is already highlighted by the 
measured value of Plt at = 0.33 GeV^ @], which is comparable to the value of Plt at 
= [0. As seen from Fig. ^ this points to an interesting structure in the region 
around 0.1 GeV^, where forthcoming data are expected from an experiment at MIT-Bates 

i- 

In the upper panel of Fig. we show the Q^-dependence of the VCS response function 
Pll - Ptt/£, which reduces at the real photon point {Q^ = 0) to the electric polarizability a. 
At non-zero Q^, Pll is directly proportional to the scalar CP a{Q^), as seen from Eq. (|66D, 
and the response function P^t of Eq- (0) contains only spin GPs. As is shown by Fig. |^, 
the nN dispersive contribution to a and to the spin GPs are smaller than the asymptotic 
contribution to a, which is evaluated for Aq, = 1 GeV. At = 0, the vrA^ dispersive and 
asymptotic contributions to a have the same sign, in contrast to (3 where both contributions 
have opposite sign and largely cancel each other in their sum. 

The response functions Plt and Pll - Ptt/£ were extracted in [^] by performing a LEX 
to VCS data, according to Eq. (^). To test the validity of such a LEX, we show in Fig. |^ 
the DR predictions for the full energy dependence of the non-Born part of the ep ejrf 
cross section in the kinematics of the MAMI experiment [Q. This energy dependence is 
compared with the LEX, which predicts a linear dependence in q' for the difference between 
the experimentally measured cross section and its BH -|- Born contribution. The result of 
a best fit to the data in the framework of the LEX is indicated by the horizontal bands 
in Fig. 1^ for the quantity [d^a — (^So-BH+Bom^ j ^^/^ where $ is a phase space factor defined 
in 0. The fivefold differential cross section d^a is differential with respect to the electron 
lah energy and lah angles and the proton cm. angles, and stands in all of the following 
for da I dkf^i^dQ'f^fjdQP ^ . It is seen from Fig. ^ that the DR results predict only a modest 
additional energy dependence up to q' ~ 0.1 GeV/c and for most of the photon angles 
involved, and therefore seems to support the LEX analysis of [Q. Only for forward angles, 
G^-™- ^ 0, which is the angular range from which the value of Plt is extracted, the DR 
calculation predicts a stronger energy dependence in the range up to q' ^ 0.1 GeV/c, as 
compared to the LEX. It will be interesting to perform a best fit of the MAMI data using 
the DR formalism, extract the two fit parameters a(Q^) and /5(Q^), and consequently the 
values of Pll — Ptt/^ and Plt respectively. Such a best fit using the DR formalism is 
planned in a future investigation. 

Increasing the energy, we show in Fig. ^ the DR predictions for photon energies in the 
A(1232)-resonance region. It is seen that the ep epy cross section rises strongly when 
crossing the pion threshold. In the dispersion relation formalism, which is based on unitarity 
and analyticity, the rise of the cross section with q' below pion threshold, due to virtual nN 
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intermediate states, is connected to the strong rise of the cross section with q' when a real 
ttN intermediate state can be produced. It is furthermore seen from Fig. ^ (lower panel) 
that the region between pion threshold and the A-resonance peak displays an enhanced 
sensitivity to the GPs through the interference with the rising Compton amplitude due to 
A-resonance excitation. For example, at q' ^ 0.2 GeV/c, the predictions for Plt in the 
lower right panel of Fig. || for = 0.4 GeV and A/3 = 0.6 GeV give a difference of about 
20 % in the non-Born squared amplitude. In contrast, the LEX prescription results in a 
relative effect for the same two values of Pit of about 10% or less. This is similar to the 
situation in RCS, where the region between pion threshold and the A-resonance position 
also provides an enhanced sensitivity to the polarizabilities and is used to extract those 
polarizabilities from data using a DR formalism. Therefore, the energy region between 
pion threshold and the A-resonance seems promising to measure VCS observables with an 
increased sensitivity to the GPs. The presented DR formalism can be used as a tool to 
extract the GPs from such data. 

When increasing the value of e, the Born and non-Born parts of the ep — epy cross section 
increase relative to the BH contribution, due to the increasing virtual photon flux factor [|I| . 
This is seen by comparing the non-Born cross section in Fig. (corresponding to e = 0.62), 
with the result for e = 0.8 at the same value of q and O^^'™', as is shown in Fig. ^ Besides 
giving rise to higher non-Born cross sections, an experiment at a higher value of e (keeping 
q flxed) also allows to disentangle the unpolarized structure functions Pll{.(]) and PTT(q) in 
Eq. (^). This will provide a nice opportunity for the MAMI-C facility where such a higher 
e value (as compared to the value e = 0.62 of the flrst VCS experiment of Ref. 0]) will be 
reachable for the same value of q. 

Recently, VCS data have also been taken at JLab [Q both below pion threshold at = 
1 GeV^ and at = 1.9 GeV^ [^, as well as in the resonance region around = 1 



GeV^ m 



The extraction of GPs from VCS data at these higher values of Q , requires an accurate 
knowledge of the nucleon electromagnetic form factors (FFs) in this region. For the proton 
magnetic FF G\j{Q'^), we use the Bosted parametrization [^, which has an accuracy of 
around 3 % in the region of 1 - 2 GeV^. The ratio of the proton electric FF to the 
magnetic FF was recently measured with high accuracy in a polarization experiment at 
JLab in the range 0.4 - 3.5 GeV^ [Q. It was found in that drops considerably 



faster with than G^j. In the region of interest here, i.e. in the 1-2 GeV^ range, the 
JLab data of Ref. pS] are well described by the flt PT| : 



0.13 {Qy + 0.028 {Qy , (69) 



where fip is the proton magnetic moment. In the following VCS calculations at = 1 GeV^, 
we use the parametrization of Eq. (^) to specify G^ (with the Bosted parametrization for 
Gm)- 

In Fig. |l^, we show the DR predictions for the ep — > epj reaction at = 1 GeV^, for 
three values of the outgoing photon energy, below pion threshold. In these kinematics, data 
have been taken at JLab and, at the time of writing this paper, preliminary results on VCS 
cross sections and GPs have been reported in |^ . For those kinematics, we show in Fig. ^ 



the differential cross sections as well as the non-Born effect relative to the BH + Born cross 
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section. It is seen from Fig. |T0| that the sensitivity to the GPs is largest where the BH + 
Born cross section becomes small, in particular in the angular region between 0° and 50". 
In Fig. |1^, we show the non-Born effect for different values of the polarizabihties. For Plli 
the calculation for the vriV dispersive contribution at = 1 GeV^ gives : 

P£^(lGeV2) = -0.3 GeV-2, (70) 

leading to the total results for Pn within the DR formalism : 

Pll{1 GeV^) = +2.3 GeV"^ for (A„ = 1 GeV) , (71) 
Pll{1 GeV^) = +4.2 GeV-^ for (A„ = 1.4 GeV) . (72) 

For PlT) the calculation for the vrA^ dispersive contribution at = 1 GeV^ gives : 

Pl^ilGeY^) = -0.9 GeV-2, (73) 

leading to the total results for Plt within the DR formalism : 

PLT(lGeV^) = -0.6 GeV-\ for (A^ = 0.6 GeV) , (74) 
PLT(lGeV^) = -0.9 GeV"^ for (A^ = 0.4 GeV) . (75) 

It will be interesting to compare the sensitivity of the cross sections to these values of the 



GPs, as displayed in Fig. |T^, to the JLab data which have been taken in this region |^ 
The deviation of the experimental values from the dispersive ttN values of (^) for Pn and 
of ( [75| ) for Plj- will provide us with interesting information, allowing to test our understand- 
ing of the electric and magnetic polarizability at this large virtuality of = 1 GeV^. 

For the same kinematics as in Fig. [1^, we compare in Fig. [TT| the DR calculation for the 
non-Born cross section with the corresponding result using the LEX. It is seen that the 
deviation of the DR results from the LEX becomes already noticeable for q' = 75 MeV, over 
most of the photon angular range. Therefore, the DR analysis seems already to be needed 
at those lower values of q' to extract GPs from the JLab data. 



In Fig. |T2|, we increase the energy through the A(1232)-resonance region, and show the 
results for the ep — >• ep'j reaction at = 1 GeV^ and at a backward angle. We display 
the calculations of the cross section and of the non-Born effect for the values in ( [7I| ) and 
(ff^) for Pll, and for the value in (|7^ for P^t- One sees a sizeable sensitivity to Pll in this 
backward angle cross section, and it therefore seems very promising to extract information 
on the electric polarizability from such anticipated data. 

Until now, we discussed only unpolarized VCS observables. An unpolarized VCS exper- 
iment gives access to only 3 combinations of the 6 GPs, as given by Eqs. (|66|-|68|). It was 
shown in Ref. |jT2| that VCS double polarization observables with polarized lepton and po- 
larized target (or recoil) nucleon, will allow us to measure three more combinations of GPs. 
Therefore a measurement of unpolarized VCS observables (at different values of e) and of 3 
double-polarization observables will give the possibility to disentangle all 6 GPs. The VCS 
double polarization observables, which are denoted by AA^(/i, z) for an electron of helicity 
h, are defined as the difference of the squared amplitudes for recoil (or target) proton spin 
orientation in the direction and opposite to the axis i {i = x, y, z) (see Ref. [|12| for details). 
In a LEX, this polarized squared amplitude yields : 
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AA^^^P = —4^ + + AMr + O(q') . (76) 

q' q' 

Analogous to the unpolarized squared amplitude (0), the threshold coefficients AAI^^I", 
AAi^_^i are known due to the LET. It was found in Refs. ]T^,P^ that the polarized squared 
amplitude AAi"^^ can be expressed in terms of three new structure functions Plj,{q), P^^(q), 
and Pjjp{q). These new structure functions are related to the spin GPs according to ||T^,[T[| : 

Pl^ = ^ Gm P''''''^' - ^ P^'''''^' , (77) 



LT 



l Q Gm P^'''''^' + Ge P^''^''^' , (78) 



2^ Qqo 



P'l^t = i^P^'''''^' - y|goP(^^'°^)^j . (79) 

While PIj, and 

P^T t)e accessed by in- plane kinematics ($ — 0°), the measurement of 
P]^ requires an out-of-plane experiment. 

In Fig. |l^, we show the dispersion results for the double polarization observables, with 
polarized electron and by measuring the recoil proton polarization either along the virtual 
photon direction (2;-direction) or parallel to the reaction plane and perpendicular to the 
virtual photon (x-direction). The double polarization asymmetries are quite large (due to 
a non- vanishing asymmetry for the BH + Born mechanism), but our DR calculations show 
only small relative effects due to the spin GPs below pion threshold. Although these ob- 
servables are tough to measure, a first test experiment is already planned at MAMI [0. 

When measuring double polarization observables above pion threshold, one can enhance 
the sensitivity to the GPs, as we remarked before for the unpolarized observables. In Fig. |T^, 
we show as an example the double polarization asymmetry in MAMI kinematics for polar- 
ized beam and recoil proton polarization measured along the virtual photon direction as 
function of the outgoing photon energy through the A(1232) region. The A(1232) reso- 
nance excitation clearly shows up as a deviation from the LEX result above about q' = 100 
MeV. 

As discussed before, VCS polarization experiments below pion threshold, require the mea- 
surement of double polarization observables to get non-zero values, because the VCS am- 
plitude is purely real below pion threshold. However, when crossing the pion threshold, the 
VCS amplitude acquires an imaginary part due to the coupling to the vrA^ channel. There- 
fore, single polarization observables become non-zero above pion threshold. A particularly 
relevant observable is the electron single spin asymmetry (SSA), which is obtained by flipping 
the electron beam helicity |[l[]. For VCS, this observable is mainly due to the interference of 
the real BH + Born amplitude with the imaginary part of the VCS amplitude. In Fig. |l^, 
the SSA is shown for two kinematics in the A(1232) region. As the SSA vanishes in-plane, 
its measurement requires an out-of-plane experiment, such as is accessible at MIT-Bates 
p6| . Our calculation shows firstly that the SSA is quite sizeable in the A(1232) region. 
Moreover, it displays only a rather weak dependence on the GPs, because the SSA is mainly 
sensitive to the imaginary part of the VCS amplitude. Therefore, it provides an excellent 
cross-check of the dispersion formalism for VCS, in particular by comparing at the same 
time the pion and photon electroproduction channels through the A region. 
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VIII. CONCLUSIONS 



In this work, we have presented a dispersion relation (DR) formahsm for VCS off a proton 
target. Such a formahsm can serve as a tool to extract generahzed polarizabihties (GPs) 
from VCS observables over a larger energy range. The way we evaluated our dispersive 
integrals using nN intermediate states, allows to apply the present formalism for VCS ob- 
servables through the A(1232)-resonance region. 

The presented DR framework, when applied at a fixed value of Q^, involves two free pa- 
rameters which can be expressed in terms of the electric and magnetic GPs, and which are 
to be extracted from a fit to VCS data. We proposed a parametrization of these two free 
parameters (asymptotic terms to a and j3) in terms of a dipole dependence, and inves- 
tigated the sensitivity of VCS observables to the corresponding dipole mass scales. 

We confronted our dispersive calculations with existing VCS data taken at MAMI below 
pion threshold. Compared to the low energy expansion (LEX) analysis which was previously 
applied to those data, we found only a modest additional energy dependence up to photon 
energies of around 100 MeV, which supports such a LEX analysis. When increasing the 
photon energy, our dispersive calculations show that the region between pion threshold and 
the A-resonance peak displays an enhanced sensitivity to the GPs. It seems therefore very 
promising to measure VCS observables in this energy region in order to extract GPs with 
an enhanced precision. 

Furthermore, we showed our DR predictions for VCS data at higher values of Q^, in the 
range = 1 - 2 GeV^, where VCS data have been taken at JLab which are presently 
under analysis. It was found for the JLab kinematics that the DR results show already a 
noticeable deviation from the LEX result even for outgoing photon energies as low as 75 
MeV. Therefore, the DR analysis seems already to be needed below pion threshold to extract 
GPs from the JLab data. We also showed predictions ai = 1 GeV^ at higher outgoing 
photon energies, through the A(1232)-resonance region, where data have also been taken at 
JLab. At backward scattering angles, we found a very sizeable sensitivity to the generalized 
electric polarizability. The two different JLab data sets, both below pion threshold and in 
the A-region, at the same value of Q"^ (in the range Q"^ — 1 - 2 GeV^) will provide a very 
interesting check on the presented DR formalism to demonstrate that a consistent value of 
the GPs can be extracted by a fit in both energy regions. 

Besides unpolarizcd VCS experiments, which give access to a combination of 3 (out of 6) 
GPs, we investigated the potential of double polarization VCS observables. Although such 
double polarization experiments with polarized beam and recoil proton polarization are quite 
challenging, they are needed to access and quantify the remaining three GPs. Using the DR 
formalism one can also analyze these observables above pion threshold. 

Finally, above pion threshold also single polarization observables arc non-zero. In partic- 
ular, the electron single spin asymmetry, using a polarized electron beam, is sizeable in the 
A-region and can provide a very valuable cross-check of the VCS dispersion calculations, as 
it is mainly sensitive to the imaginary part of the VCS amplitude, which is linked through 
unitarity to the ttN channel. 
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APPENDIX A: GAUGE-INVARIANT TENSOR BASIS FOR VCS 

1. VCS tensor basis p^^ 

In writing down a gauge-invariant tensor basis for VCS, it will be useful to introduce the 
following symmetric combinations of the four-momenta (in the notations of Sec. |I|) : 

P=\{p + p'). K=]^{q + q'). (Al) 

The 12 independent tensors p^^ entering the VCS amplitude of Eq. (H), that were intro- 
duced in [|ll|, are given by : 



pf = -(2Mz/)^^^'^ - 4g-g'P^P^ + 4Mz/ {p^'q'' + P'^g'^) , 
= -2MvQ'^g^"' - 2Muqf'q'' + 2Q'^P''q'^ + 2g-g'PV , 
pf = SP^P":^- AMv (P^Y + P'^7'^) + i AMv 75 £'^'^"'^fs:„7^ , 

= -Sg-g'P/^P^ + AMu (p^q"" + P'^g''^) + 4Mg-g' (^P'^7^ + P''7^ 
-4MV (^g''^7^ + 5^7/^) +iAMv (^q'^a'^Ka - q^a^'^K^ + q-q'a""^ 

+i4Mg-g'75£^^"^ir„7/3, 
p^'^ = (^P^g'^ - P'^g'^j^- g.g' (^P^7'^ - p-7'^j + Mz/ (^g''^7'^ - g'^7^) , 



p^'^ = Muqf'q" + - ^''^'^) - q-q' P^q^ - Mq^'q''$+ Mq-q'qf'Y 

M 



p7 



+ yQ' {q'^i" - q"Y) - ^ (g"^^""i^a - gV^°ir, + g-gV 

2Mu {P^q'' - P^g'^) - 2Mg-g' {p^Y - P'^7^) + 2MV (^g'^7^ - g^7'^ 
+i 2g-g' (^P'^a'^'^K^ + P^'a^'^'K^ - 1 2Mv (q'^'a^'^'K^ + gV^^/sT^) , 
-AMu g^'" + 2 [P^'q" + P'^g''^) + 4M g^'''$- 2M (^q'^'Y + q'' 
-2i (^g'^(T^"K„ - q'^af'^K^ + q-q'a' 



pA'- = 4 (P^g'^ + P'^g"^)^- AMu (g'^7^ + g'^7^) + 2 4g-g' 75 £'^'^"^7^,7/3 , 
p/^^ = 2Q2pMp^ + 2Mz/ P'^g'^ - 2MQ^P''Y - 2MV g''7'^ + i 2Mv q'^a'^'Ka 

+i [p^'a^^'K^ + P'^a^^'K^ - Mu a^""^ - i MQ^ 75 e'''^"^ir„7^ , (A2) 
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where we follow the conventions of Bjorken and Drell ||3^, i.e. a'^'^ = 2/2 [7^^, 7^] and in 
particular £0123 = +1- 



2. VCS invariant amplitudes Bi of Berg and Lindner 



For further reference, it also turns out to be useful to work with an alternative tensor basis 
for VCS, introduced by Berg and Lindner [|T7 . 



One starts by defining, besides the four- vectors P and K of Eq. ( |A1|) , the combination 

1 



(A3) 



and constructs from K, P, and L, the following four- vectors which are orthogonal to each 
other : 



{L-K) 



L'2 



(A4) 



One next constructs the combination of the four-vectors K and L' which is gauge-invariant 
with respect to the virtual photon four-momentum q as : 



a ■ K 

K'" = - L'^" 

q-L' 



(A5) 



which satisfies q ■ K' = 0. In terms of these four-vectors, the Lorentz- and gauge-invariant 
VCS tensor Ai'^'^ can now be written as : 



+ ^1. (b, ^75 + 56^) + \,,S. ( Br ^75 + Bslj^ 



+ 



pi2j^2 



- [B, + 5,0^) + [Bn Z75 + i?i2^ 



K^P' 



(A6) 



where Bi{i = 1, 12) are the VCS invariant amplitudes of Berg and Lindner |T7 



The invariant amplitudes defined in Eq. (11) which correspond to the tensor basis of 
Eq. ( |A2D can be expressed in terms of the invariant amplitudes Bi defined in Eq. 
These expressions read : 



Fi 



it + Q'^fP'^ 



Bi-B3 + u{B2 



B. 



+ {t + Q^P'^ 



B3 + uBi 



M 



P' 



-u%t - Q')Be + 8MuQ' 



Bo + uBw - 



t 



2(t + Q2)p'2 



5i - 53 + u{B2 - 54) 



4M 
t 



■B 



2MP'2 



12 



B. 
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1 



2M 



(t + g2)2p/2 

H Bg + vBiQ — -j^Bi2 



AM 



(i + g2)p/2 

4M 



^2 ~ B4 

' Mu 



Be , 



Mi 



(i + g2)2p/4- 

2z/2(t - Q 



2^-, 



1 + 



(t + Q2)2 



i?2 — -B4 



M 



4 



B, 



{t + Q'f 



M{t + Q2) 
32z/g2 



Br — Bs — 2 



1 + 



1 



Bn + 



8u{t-Q^) 



it + Q 



2\2 



:{-j{t' + tQ' + 2Q' 



4M2-t/ 2z/2(t-Q2^ 
2P'2 

-Bio 

4(t-g2). 

— Ba 



B. 



+ 



i?12}, 

g2(t + g2)2p'2 



g^ 4g2 



g^ 



4g2 

^4 + p72 ^5 



(2Mi/g2 



+2tBs - AQ^Bw + ^ T^B 



8M 
pi 2 



B, + ^[iAM'-t)-2P''^B, 



(t + Q 



2\2 



12 



(t + Q-'yp 



<I2 



{- 



2Mv 



+AB 



16 



10 



Bn + 



B2 - P4J + 
AM^ -t 



Mv 



-B 



12 



(t + g2)3p'2 



Bo — Bi 



8Q' 



4. 



AMu 

1 rMz/r 2z/2(t-g2) 

(t + g2)2p'2i^r i + g2 

SMu^Q^ 1 

^(i + g2)2p'2 5+p7^ 



P9 — P4 



- - U + 



2M 

2z/2(t - g^ 



{t + Q')P"B4 

tP''^-\ 2 

Pfi — — B7 



2 J 



+ 



1 + 



Aiy' 



P 



10 



i6z/g2 4Mz/g2 p g2 
Pii + 



(t + g2)2 

P2 ~ B4 

2 



t 



t + g2 4M2 



p 



12 



}. 



+ g2)2p' 
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APPENDIX B: BORN CONTRIBUTIONS TO INVARIANT AMPLITUDES 



For the invariant amplitudes Fj, defined through Eq. ([T^), one finds the following Born 
contributions , corresponding to a nucleon intermediate state in the s- and w-channel of 
the 7*p 7P process : 



F.^ 



2 

B 



F^ 



6 
^7 

fF 



^9 



7B 
10 



pB 



?B 
12 



1 ft 


M [^s — 


M2)(m- 
1 


M2) I 


M{s 


-M2)(m 


-M2) [ 


0, 


1 




2M( 


s-M2)(u-M2) 


M2(s- 


-M2)(m- 

1 


-M2) I 


4M(s - 


-M2)(m- 


- M2) V 


0, 






0, 


1 




2M(s - 


-M2)(m- 

1 


-M2) [ 

/'I 1 



kFi{Q^) + (1 + /€)F2(g')] - z/2/tF2(Q2)} 



i^i(Q') + 



8M2 



/€Fi(g2) + (i + 2/s:)F2(Q2) +zy2^^2(Q2 



(2 + «:)Fi(Q^) + F2(Q^) + -^^^F,{Q 



{s-M^){u-M^y 
1 

4M2(s-M2)(m-M2) 



-«:Fi(g^) + F2(Q^) + ^«:F2(g 



/€Fi(Q2)+F2(g2) _^,2^^^(g 



8M3(s-M2)(m-M2) 
where and F2{Q'^) are the Dirac and Pauli nucleon form factors respectively. 
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APPENDIX C: S-CHANNEL HELICITY AMPLITUDES FOR VCS 

1. Definitions and conventions 

The s-channel helicity amplitudes for virtual Compton scattering are denoted by 
Ty y and were defined in Eq. (j^). In this appendix, we express the invariant ampli- 

tudes Fi in terms of these s-channel helicity amplitudes. In addition, we quote the explicit 
results for the imaginary parts of the helicity amplitudes in the case of irN intermediate 
states. 

We work in the cm. system of the s-channel process 7*A^ —>■ jN, and all kinematical 
quantities are understood in this system. The energies of the incoming (outgoing) nucleon 
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are denoted by E {E') respectively. The incoming photon has energy go and its momentum 
q is chosen to point in the 2;-direction. The outgoing photon momentum q' is chosen to he 
in the xz plane and makes an angle 9 with the 2;-axis. We use the Lorentz gauge for the 
photon polarization vectors. For the initial (virtual) photon, the transverse and longitudinal 
polarization vectors are given by : 



.''(,.A = ±l)=(o,T-L.--!.,o), 
,.(,,A = 0)=(|,0,0,|), 



whereas for the final (real) photon, the polarization vectors are given by : 
e'^{q',X'-±l) = (0,T-Leos^,--j=,±-Lsin^). 



(CI) 



(C2) 



The initial nucleon, characterized by the momentum p and the polarization Ajv, is propa- 
gating in the negative ;2-direction. The final nucleon, with momentum p and polarization 
A^, makes an angle 180° — 9 with respect to the virtual photon, and has the azimuthal angle 
180° + (I) J* J. This leads to the following spinor conventions for the incoming and outgoing 
nucleons : 



where 



u 



(p, A^) = VE + M 



9 \ IpI v 

E+M . 



u(p, A'jv) = VE' + M 



E+M 



9V IP I v' 



(C3) 



2 



^-1 



-1 





sm ■ 



— cos : 



2 



COS ^ 



sm : 



(C4) 



2. VCS reduced helicity amplitudes 

The reduced helicity amplitudes Tj are defined by factorizing out from the helicity ampli- 
tudes T^iy .^^^ the kinematical factors in (cos f)'^^^ ' and (sin |)'^ ^ ' , with A = A — Aat 
and A' = A' — A^. The relations between the 12 independent VCS helicity amplitudes and 
the reduced helicity amplitudes {i — 1, .., 12) read : 
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Tii.^i = cos-ri, 
T,„.^,. =cos -sm-rs, 

rrs . ,e 

T_^^,. =cos-sm -rr, 

rrs ■ ^ 

rrs . ,9 

=sm-cos -rn, 



= cos-sm -r4, 

2 ' 2 Z Z 

^ 2 • ^ 2 2 

2 ^ ^ 

T 1, 1 = cos - sm -Tg, 
2' -^2 2 2 

^-i-i;oi = cos-rio, 
T^i_i;oi = cos-sm -ri2. 



(C5) 



3. Relations between the invariant amplitudes of Berg and Lindner and the VCS 

reduced helicity amplitudes 

The imaginary parts of the invariant amphtudes Fi, which enter the dispersion integrals of 
Eq. (pTSf) , are constructed from the VCS reduced helicity amplitudes Xj, which were defined 
in (|(J5D. To avoid too lengthy formulas, we display here the relations between the amplitudes 
Bi and the Tj. The relations between the Fi and the Xj are then obtained from those relations, 
and by using Eq. which expresses the Fi in terms of the Bi. 

For convenience we define the following abbreviations for kinematical factors : 



C, = l 
Cs = l + 



E + M ' 
^ y/j-M |g| 
y/s + M E + M 

E + M E' + M ' 



C2 



1 ^ , 

E + M ' 

^ v^-M |g| 
y/s + M E + M 

\q\ 

E + M E' + M ' 



With these definitions, the relations between the amplitudes Bi of Berg and Lindner and 
the reduced helicity amplitudes Xj are given by : 



ViE + M){E' + M)_^^_^^ \q\-q,cos9 



C2 



e ( 



Ti + COS 2 ^2 - sm - \Ti + T-j 



9 o 9 



Tg - sin^ ^ ^6 + cos^ ^[n + Ts] + 



t + Q^ 



■ 2^ ( 2 ^ 

-I a^^^ - Tg-COS -Tn 

g I — go cos 9 z\ 2 



Q ( 



9 



(C6) 



I -I ^cos-- rio-sm--ri2 

|g I — go cos 9 z\ 2 



Bo 



^{E + M){E' + M) |g| -goCOS^ 



4v^|g| 



t + Q' 



(C7) 
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+ C4 



2 9 2^/ \ 2V2Q 

n + cos -r,-sin 3 + " 1,-1 - cos ^ 



sm - [tq - cos - Til 



2- , 2^/ ^ A , 2x/2<9 



Ts — Sin^ - Tq + COS 

2 



Or 



9 



q \ — go cos 9 



cos - Tio - sm - T12 



^{E + M){E' + M) 



C2 



8^|g||g' 
9 



9 



2^ ■ 2^ ( \ 
Ti + cos -Ta + sm ^[ji + ^T) 



r5 + sin^ ^ '^6 + cos^ ^ ( T3 + Ts 



(C8) 



V(^ + M)(E' + M) 



S^s\q\\q'\ 



C3 
+ C4 



2 



2" ■ 2^ ( \ 

Ti + cos - r2 + sm - (^T4 + j 



r5 + sin^ ^ '^6 + cos^ + n 



(C9) 



V(^ + M)(^' + M) s |g |2 |gf sin^ ^ 



4i 



C, 



(t + g2)2 
+ 



go) (t5 - sin^ ^ re) + (|g I + go) sin^ | (^Tg - Tg) 



• 2^ 
Tio + sm -T12 



) 



, 9 \ 9 1 \ 

(1^1 + ?o) (ri - cos^ - Taj - (|g | - go) cos^ - [r^ - Tyj 

-V^Q^Tg + COS^ ^Tii 



(CIO) 



V(E + M)(E' + M) |g'P sin^^ 



2^2 



C5 sin^ - 



(t + Q 

/ 9 \ 9 / \ 

+ go) (ti - cos^ i^r2j -{\q\- go) cos^ - (^r4 - rrj 

9 



Tg + COS - Til 



(Cll) 
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+ Cq cos^ 



/ 9 \ 9 / 

\ - Qo) - sin^ - T6 j +{\q\+ Qo) sin^ - [r^ - rs 

9 



Tio + Sm -Ti2 



^y{E + M){E' + M) s\q\'^\q'\^sm^9 



C5 sin^ 



9 / o9 



— Ce cos^ - 



{\q \ + go) sin^ - (t5 - sin^ - re) + (|g | - go) cos^ 2 (^^ ~ 

(5 COS^ ^ (^Tio + Sin^ ^ Ti2^ 

{\q \ - qo) cos^ ^ (^n - cos^ ^ - (|g | + go) sin^ ^ (|t4 - T7 

( 



—V2Q sin^ 2 ( ^9 ^ 



V(£^ + M)(£;' + M) Vs|g| Igfsin^^ 



(|g I - go) cos^ ^ (^n - cos^ ^ ra) - (|g | + go) sin^ ^ (^r4 - ry) 



-V2Q sin^ ^ ^Tg + cos^ | rn 
(|g I + go) sin^ ^ (^rs - sin^ ^ re) + (|g | - go) cos^ ^ (^ts - Tgj 

-V2Q COS^ ^ ^Tio + Sin^ ^ Ti2^ 



+ C6 



y(^ + M)(^' + M) 



8 (t + g2) 



C2 



Ti + cos^ 2^^~ 2^^^ + '^^ 



(V^-M) Ig'lsin^^ 

^ / \ 1^1 — Q'o cos ^ / 



V2Q COS 



2^ 

Tg - COS - Til 



2 y'^ — 2 



• 2 ^ 2 

T5 — sm - Te + cos 



+ 



go cos 9 



\pl Q sin 



2 e 



Tio - sm - T12 



s{t + g2) 



I g I sin 9 
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+ C4 



Ti + COS - T2 



Sm- + 



\q \ — go cos 9 / o 

Tg - COS - Til 



Q COS 



2£ V" 2 



2 ^ 2 ^ 

Ts - sm 2 ^6 + COS - 



q\ — go cos 6* / 9 ^ 

^3 + ^8 ) + ' . ^ ^ (^rio - sm - ri2 



y2Q 



sm 



[-e')Bn 



y/{E + M){E' + M) 



y/s |g I |g'p sin^ 6^ 



C5 sin^ - 



Cfi cos^ - 



rg - sin^ 2 2 ^ ~ ^® 

Ti - cos^ ^ ^2 - sin^ ^ ( 7-4 - rr 



(C16) 



|g I — go cos6' / 



V2Q 



sm 



^(^Tio + sm -ri2 



|g I — go cos 9 
V2Q 



cos 



2 e 



Tg + COS - Tii^ 



;-e^)5i2 



C5 

+ C6 



8t 

Ti - cos^ 2 ~" ^^^^ 2 ' ~ ^'^ 



Ts - sin^ ^ ^6 + cos^ ^ ( ^3 - T-8 ) - 



(C17) 



|g I — go cos6' 



Q cos 



„ , To + cos - Til 

2 ^ V 2 



g I — go cos6' / 2 ^ \ 

f-^.^r. ^10 + sm - ri2 
V2Qsin2^ V 2 /_ 



Eqs. ( [C6| - |C17| ) together with the equations in ( [AT] ) eventually allow to express the in 
terms of the VCS helicity amplitudes. 



4. Unitarity relations between the VCS reduced helicity amplitudes and pion photo- 

and electroproduction multipoles 

If we write down the unitarity equations for the VCS helicity amplitudes and consider 
only vrA^ intermediate states, then the imaginary parts of the VCS helicity amplitudes can 
be expressed in terms of the 7*A^ — > vrA^ times — > vrA^ multipoles : 

Imn = -Svrg^v^ ^^(2/ + 2) {Ai^[(f) A\^{^) + A(i+i)_(Q2)A^,+i)_(0)] 

xF(-/;Z + 2;l;sin2^), 
Imr2 = -87rg.v^5^^^^±M±^ [i?,+ (Q2)i?;+(0) + i?(,+i)_(g2)i?(V,)„(0)] 

xF(-/ + l;/ + 3;l;sin2-), 
Imrg = -87rg.v^5^/(/ + 1)(/ + 2) {Ai^{(f)B\^{^) + A(,+i)_(Q2)5(Vy_(0)] 

Q 

xF(-/ + l;/ + 3;2;sin2-), 
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xF(-/ + l;/ + 3;3;sin2-) 



Q 

xF(-/;/ + 2;2;sin2-; 



/^(/ + l)^(/ + 2)- 
12 



[i?,+(g^)i?r+(o)-i?(,+i)_(g2)i?(l 



('+!)- 



a 

xF(-/ + l;/ + 3;4;sin2-; 



Imr7 = — 87rg7r-\/s 



/(/ + l)2(/ + 2) 



[A,4g2)i?r4o) - A(,+i)_(g2)i?(Vi)_(o)] 



xF(-/ + l;/ + 3;3;sin2-), 
Imrs = -87rg^v^^/(/ + + 2) [Bi^{Q^)A*^{Q) + i?(,+i)4g2)A^,^y_(0)] 

xF(-/ + l;/ + 3;2;sin2-), 



Imrg = STrg^v^^^ + 1)^ [U+{Q^) A*i^{Q) + L(,+i)_(g2)A^^ 



9o 



xF(-Z;Z + 2;2;sin2-), 
Imrio = STrg.v^^^ + 1)^ [L,+ (g2)A*+(0) - L(,+i)4g2)A^^^^)_(0)] 

^0 , 



xF(-/;/ + 2;l;sin2^), 



?0 



xF(-/ + 1;/ + 3; 2; sin 



Imri2 = Mq.j,^/s— 



2" 



[L,+(g2)sr+(o) + L(,+,)_(g2)i?(Vi)4o)] 



xF(-/ + l;/ + 3;3;sin2-), 



(C18) 



where are the reduced hehcity amphtudes defined in Eq. (pSf). In Eqs. (|C18D , is 
the pion cm. momentum in the intermediate state, and F is a hypergeometric polynomial 
defined as : 
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^, , , ^ ab X a(a + l)b(b + 1) x'^ 

F{a; b;c;x) = l + -- + ^ \ + .. 

c 1! c(c + 1) 2! 



(C19) 



In Eq. ( |(J18| ), the transverse multipoles Ai±, Bi±, and the longitudinal multipoles L/± are 



defined as in Ref. 38 
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FIGURES 




FIG. 1. (a) FVCS process, (b) BH process. 
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FIG. 2. The Mandelstam plane for virtual Compton scattering at = 0.33 GeV'^. The 
boundaries of the physical s-channel region are Q-y^ = 0° and B^^ = 180° for > 0, the n-channel 
region is obtained by crossing, v — > — i/. The curves for = 0° and O^^ = 180° intersect at 
u = 0, t = —Q^, which is the point where the generalized polarizabilities are defined. 
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FIG. 3. Results for the spin- flip GPs excluding the 7r''-pole contribution in different model 
calculations, as function of the squared momentum transfer. The full curves correspond to the 
dispersive vrA^ contribution. The dashed curves show the results of 0{p^) HBChPT (H, the 
dashed-dotted curves correspond to the predictions of the linear u-model , and the dotted curves 
are the results of the nonrelativistic constituent quark model ||2^ . Note that the constituent quark 
model (CQM) results for p(oi:Oi)i and p(ii.02)i ^^.^ multiplied (for visibility) by a factor 100. 
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FIG. 4. Results for the spin-flip GPs as function of the squared momentum transfer. The 
dashed curves correspond to the dispersive ttN contribution, the dotted curves show the Tr'^-pole 
contribution, and the full curves are the sum of the dispersive and Tr'^-pole contributions. For 
comparison, we also show the 7r°-pole contribution when setting the Tr^^*j form factor equal to 1 
(dashed-dotted curves). Note that p(oi'Oi)i has no TT -pole contribution. 
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( GeV^ ) 

FIG. 5. Theoretical estimates of the asymptotic contribution Ff^ : DR calculation Q of the 
7*7 vrvr — > NN process as described in the text in Sec. VI B| (full curve); linear a-model 
(LSM) calculation with nio- = 0.5 GeV (dotted curve) and nio- = 0.7 GeV (dashed-dotted 
curve). The dashed curves are dipole parametrizations according to Eq. (p3[), which are fixed to 
the phenomenological value at = and are shown for two values of the mass-scale, = 0.4 
GeV (upper dashed curve, nearly coinciding with full curve) and Afj = 0.6 GeV (lower dashed 
curve) . 
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FIG. 6. Results for the unpolarized structure functions Pll — Ptt/^ (upper panels), and Plt 
(lower panels), for e = 0.62. In the upper left panel, a comparison is shown between the dispersive 
vrA^ contribution of the GP a (full curve), the dispersive vrA^ contribution of the spin- flip GPs 
(dashed curve), and the asymptotic contribution of a according to Eq. ( |62| ) with = 1 GeV 
(dotted curve). The upper right panel displays the total result for Pil — Ptt/^ (sum of the three 
contributions on the upper left panel) for Aq, = 1 GeV (full curve) and A^ = 1.4 GeV (dashed 
curve). In the lower left panel, a comparison is shown between the dispersive vrA^ contribution of 
the GP f3 (full curve), the contribution of the spin- flip GPs (dashed curve), and the asymptotic 
contribution of /3 according to Eq. ( ^3| ) with A^ = 0.6 GeV (dotted curve). The lower right panel 
displays the total result for Plt, for A^ = 0.7 GeV (dotted curve), A/3 = 0.6 GeV (full curve), 
and Ap = 0.4 GeV (dashed curve). The RCS data are from Ref. [^, and the VCS data at 
Q2 = 0.33 GeV^ from Ref. ||. 
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FIG. 7. {d^a — d^a^^~^^°'^'^)/^q' for the ep epj reaction as function of the outgoing-photon 
energy q' in MAMI kinematics : e = 0.62, q = 0.6 GeV, and for different photon cm. angles Q^I^'. 
The data and the shaded bands, representing the best fit to the data within the LEX formaUsm, 
are from Ref. [Q]. The sohd curves are the DR results taking into account the full q' dependence 
of the non-Born contribution to the cross section. The asymptotic contributions are calculated 
according to Eqs. (|53| , [6^ ), with = 0.6 GeV and A„ = 1 GeV, respectively. 
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FIG. 8. Upper panel: The differential cross section for the reaction ep — > ep7 as function of 
the outgoing-photon energy q' in MAMI kinematics : e = 0.62, q = 0.6 GeV, and for 0!^'™' = 0°, 
in plane = 0°). The BH + Born contribution is given by the dashed-dotted curve. The total 
DR results are obtained with the asymptotic parts of Eqs. (53, 62), using a fixed value of Aq, = 1 
GeV and for the three values of A^ as displayed in the lower right plot of Fig. ^ i.e. A^ = 0.7 
GeV (dotted curve), A^ = 0.6 GeV (solid curve), and A/? = 0.4 GeV (dashed curve). Lower 
panel: Results for (S'a — ci^(T^^^^°'''^)/<I>q' as function of q' . The DR calculation taking into 
account the full energy dependence of the non-Born contribution (thick curves) are compared to 
the corresponding results within the LEX formalism (thin horizontal curves). The curves in the 
lower panel correspond to the same values of A^ and A^ as in the upper panel. The data are from 
Ref. a. 
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FIG. 9. Same as Fig. | but for e = 0.8. 
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FIG. 10. Left panels : The differential cross section for the reaction ep — > ep-y as function 
of the photon scattering angle and at different values of the outgoing-photon energy q' in JLab 
kinematics. Right panels : ratio of cross sections {da — da^^~^^°^"')/da^^~^^°^"' . Dashed-dotted 
curves on the left panels : BH+Born contribution. The DR results are displayed (on both left and 



right panels) with the asymptotic terms parametrized as in Eqs. (62, p3), using the values : A, 



1 GeV and A/s = 0.6 GeV (full curves), Aa = 1 GeV and Ap = 0.4 GeV (dashed curves), Aq 
1.4 GeV and A^ = 0.6 GeV (dotted curves). 
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FIG. 11. (d^o" — d^cj^^"''^°™)/<5q' for the ep — > epj reaction as function of the scattering angle 
and at different values of the outgoing-photon energy q' in JLab kinematics. The solid curves 
correspond to the DR calculation with the full energy dependence of the non-Born contribution 
to the cross section. The dashed-dotted curves are the corresponding results obtained from the 
LEX. The asymptotic contributions have been calculated with the parametrizations in Eqs. (62, 
53D, using A„ = 1.4 GeV and = 0.6 GeV . 
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FIG. 12. Upper panel: The differential cross sections for the ep epj reaction as function 
of the cm. energy W in JLab kinematics : e = 0.95, = 1 GeV^, and for fixed scattering angle 
9^:^- = —160°, in plane = 0°). The BH + Born contribution is given by the dashed-dotted 
curve. The total result including the non-Born contribution is shown for = 0.6 GeV and for the 
two values : Aq, = 1 GeV (full curve), and Aq, = 1.4 GeV (dashed curve). Lower panel: results for 
[d^a — d^c7^^"'"^°™)/$q' as function of W. The curves in the lower panel correspond to the same 
values of A^ and A^ as in the upper panel. 
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FIG. 13. VCS double-polarization asymmetry (polarized electron, recoil proton polarization 
along either the z- or x- directions in the cm. frame) in MAMI kinematics as function of the photon 
scattering angle. The dotted curves correspond to the BH+Born contribution. The dispersion 
results for the total BH+Born+non-Born cross section are shown for the values of the mass-scale 
Aa = 1 GeV, = 0.6 GeV (full curves) and A^, = 1 GeV, A^ = 0.4 GeV (dashed curves). To see 
the effect of the 7r°-pole contribution, we also show the results for the values A^ = 1 GeV, A^ = 
0.6 GeV, when turning off the 7r°-pole contribution (dashed-dotted curves). 
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FIG. 14. Upper panel : VCS double-polarization asymmetry (polarized electron, recoil pro- 
ton polarization along the z-direction in the cm. frame) in MAMI kinematics (same value of q 
and e as in Fig. |l^) as function of the outgoing-photon energy at a fixed photon scattering angle 
Q'^J^' = —50°, in plane {(j) = 0°). The middle panel is the corresponding difference of polarized 
cross sections and the lower panel is the non-Born contribution to the corresponding polarized 
squared matrix element (according to Eq. (|76|)). The dotted curves correspond to the BH+Born 
contribution. The dispersion results for the total BH+Born-l-non-Born cross section (full curves) 
are calculated using the values = 1 GeV and A^ = 0.6 GeV. The dashed curves are the corre- 
sponding results obtained from the LEX. To see the effect of the vr'^-pole contribution, we also show 
the results of the dispersion calculation, when turning off the vr'^-pole contribution (dashed-dotted 
curves) . 
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FIG. 15. Electron single spin asymmetry (SSA) for VCS at = 0.12 GeV^, for two kinematics 
in the A(1232) region : W = 1.232 GeV, e = 0.75 (upper plots) and W = 1.17 GeV, e = 0.81 
(lower plots). In both cases the SSA is shown as function of the photon scattering angle for two 
out-of-plane angles as accessible at MIT-Bates |^^. The full dispersion results are shown for 
the values : Aq, = 1 GeV, A/3 = 0.6 GeV (fuh curves), A^ = 1 GeV, A/? = 0.4 GeV (dashed curves), 
and Aa = 1.4 GeV, A/3 = 0.6 GeV (dotted curves). 
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